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Abstract 

The expression for additional subgap current in the presence of electron- 
phonon interaction is derived. We show that the phonon assisted tunneling 
leads to the appearance of peaks on current-voltage characteristics at the 
Josephson frequencies corresponding to the Raman-active phonons. The re- 
lation of the obtained results to the experimental observation are discussed. 
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It is well known that the superconducting Josephson current at nonzero voltage V across 
the junction oscillates at a so called Josephson frequency uj [|l|] : 



hcu = 2eV (1) 

were V is the voltage on the contact. 

This alternating (AC) supercurrent can generate time dependent longitudinal or trans- 
verse electromagnetic fields within the junction. Nonlinear character of this phenomenon 
leads to the appearance of many peculiarities in the current-voltage characteristic (CVC) 
of a Josephson contact 0. As an example, the CVC can exhibit current peaks, or multiple- 
valued CVC with different branches can appear at some voltages. These effects result from 
the coupling of the Josephson oscillation to electromagnetic cavity modes of the junction 
structure Hf| . Voltages are connected with the corresponding mode frequencies U{ by the 
Josephson relation Eq. (|IJ). The possibility of the generation of phonons in Josephson junc- 
tions have been predicted by Ivanchenko and Medvedev ||. The authors have considered 
the excitation of the long wavelength acoustic resonance modes in the dielectric layer of 
the contact, which can influence the shape of the CVC of the junction in the same way as 
excitation of electromagnetic cavity modes. To our knowledge, the resonant current peaks in 
the CVC and the branching of the CVC, corresponding to the acoustic resonances predicted 
in ||, have never been observed experimentally. The mere fact of the phonon generation 
in classical Josephson junctions, caused by the quasiparticle current or by the AC Joseph- 
son current, was repeatedly found experimentally H-f|. This generation was observed by 
measuring phonon emission using phonon spectrometer. Two mechanisms of this generation 
were discussed: the direct conversion of electromagnetic waves into phonons and phonon- 
assisted tunneling processes || . Which of the mechanisms better fits experimental situation 
have not been established at full length. 

Recently specific subgap structures in the form of current peaks (resonances) have 
been observed flfLOfl in the CVC's of intrinsic Josephson junctions in B^S^CaiC^Og and 
Tl 2 Ba2Ca 2 Cu30io with the tunneling current in the c-direction, perpendicular to the copper- 
oxygen layers. The obtained results have been explained ||lT| , |T2"| by the coupling of the optical 



phonons in the dielectric layers being the ionic crystal with the AC Josephson current (or, 
more precisely, with the longitudinal electric field existing inside the dielectric layer). The 
displacement current contribution (D) to the total current was described in the framework 
of a simple macroscopic dielectric function for ionic crystals e(u). It was shown that peaks 
in the CVC appear when the voltage satisfies the condition 

2eV = hu h0 (2) 

where c^lo is the frequency of the longitudinal optical (LO) phonons with the wavevec- 
tor q = 0. These results can be easily obtained from the equations of Refs. ||!§ for the 
electromagnetic waves generation if one substitutes the displacement current D in the form: 

D(t) = i fduojs(oj)E(oj)e iu4 (3) 

The peak appearance corresponds to zero value of the dielectric function, because s{ujlo) — 
0. 
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A fine structure in the CVC of a single Bi 2 Sr 2 CaiCu208 break junction have been re- 



ported earlier [13,14] which could have some resemblance with the subgap structure in 
BSCCO stacks P-pll- The results of the detailed investigations of this structure are pre- 



sented in the preceding article [fLq| . The essential result of this work is the observation of 
the local branching of the CVC and the dips in the dl/dV characteristics at the voltages V 
satisfying the condition 

2eV = huj K (4) 

where hu-n are the energies of the Raman-active phonons. 

The explanation for the observed subgap structure given in Refs. [|l"l~| , |T2|| is based on direct 
conversion of the AC longitudinal electric field, existing inside the junction into LO phonons 
carrying electric field also inside the junction. But the direct conversion of longitudinal 
electric field into Raman-active phonons is forbidden due to selection rules. Moreover, the 



very existence of the well defined dielectric layers in the breakj unctions studied in Ref. (15 
is very improbable. Thus "electromagnetic" explanation can not be applied to the case 
of Raman active phonons, but such phonons can be emitted due to the electron-phonon 
interaction in the banks of the junction and can interfere with AC Josephson current. 

Here we apply the ideas developed in the work |J to explain the results obtained in Refs. 
13HT5|. For this purpose we generalize the approach used in || to include interaction of 
electrons with all phonons and not only with acoustical ones. This generalization is very 
important for the high-Tc materials, because the range of Josephson frequencies exceeds the 
range of the main part of the phonon spectra. The total Hamiltonian can be written for this 
case as 



/ 



"Htotai = 7ii + + ( X! Tk P a>i a bpa + h.c. 



T k P a L \o U q x {CgX + C± qX ). + h.C. 



k ,p,<J 
q, A 



(5) 



Here the terms "Hi and Ti,2 describe electrons on the two sides of the tunneling junction with 
the operators and b pcr correspondingly. These Hamiltonians include an interelectron 
interaction leading to superconductivity. We do not consider here this interaction explicitly. 
It can be the electron-phonon interaction or any other. The only important assumption 
made is that this interaction leads to the appearance of superconductivity and the Josephson 
current exists between the two superconductors. We will use the BCS formulas for Josephson 
current for the case of isotropic order parameter in the main part of the paper. Qualitatively, 
results are not very sensitive to the symmetry of pairing. It is supposed also that these 
Hamiltonians include the interaction of the electrons with the bias voltage 



k.cr 



and -eV 2 {t)J2 b P> 



where 



V l {t)-V 2 {t)=V{t) 



(6) 



(7) 



and V(t) is the bias voltage, applied to the junction. Since we are interested in direct 
current changes, we can use the usual approximation simplifying the problem. Namely, we 
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neglect electrodynamic selfinteraction in the junction, which is quite reasonable for weak 
AC component of the current (Josephson current) 0. Then the bias voltage is equal to a 
given constant external bias V(t) = V. The third term in Eq. (|5]) is the usual tunneling 
Hamiltonian and Tf. p is a matrix element for the tunneling of an electron. The fourth term 
in Eq. (|^) describes the phonon assisted tunneling of electrons. The value a q \ stands for 
the electron-phonon coupling and the operators c q \ and ct qX describe phonons with the 
given momentum q and the polarization A. A similar expression for this term was used 



earlier by Kleinman [[Uj to calculate the contribution of the phonon assisted tunneling to 
the quasiparticle current between superconductors. We would like to emphasize here that 
the existence of such a term in the total Hamiltonian does not depend on the origin of the 



superconductivity in HTSC systems. There is some experimental evidence |L7j that the 
HTSC material possesses considerably strong electron-phonon interaction even though it 
does not govern totally the superconducting properties of HTSC systems. 
Introducing the expression for the total current as 

/ = e(Aq) (8) 
where N\ is the number of electrons in one side of the junction and e is the electron charge 

^i=E fl ^ (9) 

fc,(T 

We can write for the current the following expression (h — 1) 

I=-ie{[N u H]). (10) 

The symbol (■ • •) means the averaging over the nonequilibrium statistical ensemble. The 
last two terms in Eq. (^J) give the nonzero contribution to this commutator. Substituting 
these terms into Eq. (|10D we get 



/ = 2elm W T pk (a+ a (t)b pa (t)) + £ T pk a qX (at(t)b pa (t) (c qX (t) + c± qX (t))) . (11) 
After long but simple calculations following close to that done in the work |5|] we get 



T (t) = y- + X> 9 a ( c «aW + ct qX (t)) J [I s (t) + In] (12) 

where Is(t) describes the contribution of the superconducting Josephson current (AC part) 
to the total current and In (DC part) gives the contribution of the quasiparticle current. 
The Josephson current can be written in the form 

I s (t) = Re { I v } sin 0(t) + Im \l p } cos <j>{t) (13) 



where 

r 

= - 

h 



^ - ^ (M) 
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The value <j>(t) is the order parameter phase difference between the two sides of the junction: 



(f)(t) = Ujt. (15) 

We shall not reproduce here the expressions for the values I p (%f) and In which can be 
found in the works P,f4],|18|l . It is only important for us that Im {I p {ujj)} and In(V) are very- 
small at low temperature T and at small voltage V. For the BCS model Re {I p (u)} at small 
voltages has the form 



Re 



u u 4 



, x ttA(T) , A(T) 
1(0) « M tanh ■ 



2 eifr 



2T 



(16) 



Here A is the superconducting gap, e is the electron charge and is the resistivity of the 
contact in the normal state. As V increases the current amplitude Re{/ p (o;)} increases also 
and at V = 2-A it has the Riedel logarithmic singularity |0| (qualitatively originated from 
the density of states singularity at the gap edges). 

To calculate the total current in the presence of the phonon assisted tunneling we should 
calculate the value (c qX + cj_ A ). Using the Hamiltonian @ and adding the free phonon part 



q,X 



(17) 



we can obtain the equation of motion in Heisenberg representation for the phonon operator 
( c q\(t) + c - q \(t))- F° r the averaged values this gives: 



^x(t) (c qX (t) + ct qX (t)) = -2Re J2 for a qX T kp ((a+ a (t)b pa (t) 



k,p,a ' 



(18) 



Here T> qX (t) is operator inverse to the phonon Green function 



V~x\t) = 



d 2 



d 



(19) 



where uo q \ are the phonon frequencies and 7 is the small phonon relaxation rate. The right- 
hand side of Eq. ([TJ) appears due to phonon-assisted tunneling processes and plays the role 
of an external source for lattice vibrations. (For the equilibrium case (c q \{t) + c^ qX (t)) = 
). Following again the work || we can calculate the contribution to this function from 
oscillating Josephson current: 



£ jdr ^ r a qX T kp (at(t)b pa (t) 

k,p 



"OA 

e 



(20) 



Since we have neglected the space dependence of the phase <p(t) only the q = component 
gives a nonzero contribution in Eq. (|20|). 

The solution of Eq. (pT8|) can be easily found now: 



' 2lUlX RehoxMl P (^-)^' 



(21) 



5 



with phonon Green function V \ (u): 



VoxH = - ^ X , (22) 

w - W 0A - W l 

Now we consider for simplicity the case of low temperature T, so only the terms with Re 



{I p } should be taken into account. In this case Eq. (21) can be written as 



cox{t) + Cqa(*)) = 2 a ° X ^ 0X Re/p (^j^j Im {V ox (uj)}shiUjt - Re{£> A (uj)} cosujt (23) 

Substituting Eq. (|23|) into Eq. (|12] ) and averaging the oscillating terms over the time, we get 
the expression for the excess current AI(V). Using the explicit form of U \ (cjj) (Eq. (|2"2])) 
we finally obtain : 

AI{V) = V] a2 ° xUJ ° X ^{ Re/ p(^/ 2 )} 2 (24) 
a e (wJ-wg A ) 

where cuj is defined by Eq. (JTJ1). This expression shows that there are peaks in the CVC at 
the voltages corresponding to the condition 

2eV = hujQ X (25) 

where ujq\ are the frequencies of the optical phonons with q = 0. The value d/(V)/dV has 
correspondingly the dips very close to voltages satisfying Eq. (p5|). The difference between 
the peaks positions and the dips is of the order of a small value 7. Note, that contrary to 
Eq. (HD, now Eq. fl25|) includes all optical phonons (not only LO) and these are the phonons 
in superconductors themselves and not in the intermediate dielectric layer. The value of 
the electron-phonon coupling a^ x determines how great the contribution from each of the 
phonon branches is. This value is equal to zero for some phonon modes with q = due to 
the symmetry reason. This is the case, for example, for infrared-active modes. As it is well 



known [^0| the coupling constant of electrons with the Raman-active phonons is nonequal 



to zero even for q = 0. Moreover, as it was shown by the first principle calculations |2T[ 



the constant of coupling with Raman-active phonons is considerably large. This is why the 
peaks mainly connected with such phonons are observed in the experiments made on break 



junctions [15 



The current excess can be find for all temperatures T and voltages V numerically using 
Eq. (|12D and Eq. (|2l|) obtained in this work. In the present work we would like only to 
mention some clear qualitative consequences of such calculations. The peak amplitudes - 
or the dips in the dI(V)/dV characteristics — will demonstrate the Riedel peculiarities: 
the amplitude will increase for the values of eV/h close to 2 A. This subgap structure 
will disappear for eV/h, larger than 2 A, because the value I P {V) quickly goes to zero for 
eV/h > 2 A. The amplitude of the peaks decreases with decreasing of the critical Josephson 
current as (0). We would like to emphasize here that the amplitude of the peaks predicted 
in the work (TT]|T2|] is also proportional to 1% in contradiction with the statement announced 
in these works. These amplitudes will also change in the external magnetic field due to the 
change of the value I c . 
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In conclusion, simple microscopic theory describing the phonon assisted Josephson tun- 
neling was derived. It was shown that there exists subgap excess current. The current excess 
demonstrates peaks at voltages V such as 2eV/h is equal to the frequencies of the Raman- 
active phonon modes. The amplitudes of the peaks are proportional to If, the square of the 
Josephson critical current. All these results are in good agreement with the experimental 
work pf. 

The case of anisotropic pairing (e.g. ci-pairing) will not change essentially the obtained 
results, till the AC Josephson current exists. The functions I p (co) and Is(w) will be cer- 
tainly different from that of given in the Ref. P,P||r5|] , but all other equations will not be 
changed. The only necessary things are the existence of the AC Josephson current and the 
electron-phonon interaction. Therefore the position of the subgap structures on the CVC 
is determined in any case by the frequencies of Raman-active modes. Of course, the form 
of this structure and its temperature dependence can differ from those for isotropic pairing. 
The problem of numerical calculations of CVC for T 7^ in different pairing models is now 
under investigation. 
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